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Many-body effects resulting from strong electron-electron and electron-phonon interactions play a 
significant role in graphene physics. We report on their manifestation in low B field magneto-phonon 
resonances in high quality exfoliated single-layer and bilayer graphene encapsulated in hexagonal 
boron nitride. These resonances allow us to extract characteristic effective Fermi velocities, as high 
as 1.20 X 10® m/s, for the observed “dressed” Landau level transitions, as well as the broadening of 
the resonances, which increases with Landau level index. 


Scanning confocal Raman microscopy has emerged as a 
key tool for studying the unique properties of graphene. 
In recent years, Raman spectroscopy has prov en to be 
highly useful not only to identify graphene^, but also 
to extract information on l ocal dopin^^H^, straiiP^ini^ 
and lattice temperatnr d^^ l ^^ l Even more insights can 
be gained when combining Raman spectroscopy with 
magnetic fields. In perpendicular B fields, electronic 
states in graphene condense into Landau levels (LLs) 
which can interact with lattice vibrations. Changing 
the magnetic field allows for both (i) tuning the Lan¬ 
dau damping of the highest optical phonon mode at 
the r point (G mode) and (ii) shifting the Raman line 
when the G mode is resonantly coupled to energeti¬ 
cally matched LL transitiond^^^^. These coupled modes 
are known as magneto-phonon resonances (MPRs) and 
they provide a unique way to study electron-phonon 
interaction and many-body physics. Although differ¬ 
ent MPRs have been measured in magneto-Raman ex¬ 
periments on graphene on graphite^^^^ and multilayer 
graphene on SiCP^ and SiOj^, for exfoliated graphene 
only a single MPR at around B = 25 T has been ob- 
serve(PiH23j. While direct electronic LL-excitations for 
single- to penta-layer graphene have been recently stud¬ 
ied for suspended graphen^^, revealing interesting in¬ 
sights into Landau level physics, the interaction of Lan¬ 
dau levels with the G mode at low magnetic fields for 
exfoliated single-layer and bilayer graphene has not yet 
been observed and investigated. 

Here we show magneto-Raman measurements on four 
different systems: (i) single-layer graphene (SLG) on 
Si02 covered with hexagonal boron nitride (hBN), (ii) 
SLG encapsulated in hBN, (iii) electrically contacted 
SLG encapsulated in hBN, and (iv) bilayer graphene 
(BLG) encapsulated in hBN. The high quality of our 
hBN-graphene-hBN sandwich devices allows us to ob¬ 
serve magneto-phonon resonances down to 2.1 T for SLG. 
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FIG. 1. (a) Schematic illustration of a graphene-hBN het¬ 

erostructure on Si02 and on hBN. (b) and (c) Optical images 
of fabricated graphene, i.e. SLG (b) and BLG (c) samples. 
The scale bars are 10 pm. The stars indicate where the Ra¬ 
man spectra in panels (d), (e) and (f) have been taken, (d) 
and (e) Characteristic Raman spectra of graphene on hBN (d) 
and on Si02 (e). (f) Characteristic Raman spectrum taken on 
the sample depicted in panel (c). The inset shows a close-up 
of the 2D peak. The colored curves represent Lorentzian fits 
to the four sub-peaks with the yellow line representing the 
sum of the four fitted peaks. 


We are able to quantitatively study the influence of the 
electron-phonon interaction on the G mode at low mag¬ 
netic fields by pinning the charge carrier density in our 
samples. In particular, we extract characteristic lifetimes 
for the observed Landau level excitations as well as high 
values of the effective Fermi velocities, a hallmark of 
electron-electron interaction effects, which appear to be 
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LL transition and/or B field dependent. Finally, we re¬ 
port on MPRs in bilayer graphene. Here the extracted 
effective Fermi velocities are lower, showing that many- 
body effects are less pronounced as compared to single¬ 
layer graphene. 

A schematic illustration of a typical sample is shown in 
Figure [^, where a graphene flake is partly deposited on 
Si02 and partly on hBN. The upper side of the graphene 
flake is completely covered with a second hBN flake. This 
type of sample gives us invaluable capability to com¬ 
pare the material parameters of the same graphene flake, 
where one surface is in direct contact with two kinds of 
substrates. We employ a dry and resist-free fabrication 
method similar to Ref.|26l where we pick-up an exfoliated 
graphene flake with an hBN flake and deposit it onto the 
hBN-Si02 transition area of the substrate. This pro¬ 
cedure has been shown to produce hi gh-qu ality devices, 
as proven by transport measurements^^^. An optical 
image of such a structure with single-layer graphene is 
shown in Figure [^. In Figures and[^ we show two 
Raman spectra from the different substrate regions (see 
labeled regions and stars in Figure [^). In both cases the 
characteristic hBN, G and 2D peaks are observed. The 
single Lorentzian shape of the 2D line is characteristic for 
SLG. The 2D line also contains information on the un¬ 
derlying substrate. For graphene encapsulated in hBN, 
we regularly find values of the full width at half maxi¬ 
mum (FWHM) down to 16 cm“^, while Si02 supported 
graphene shows values above 22 cm“^. An optical image 
of a bilayer graphene flake encapsulated in hBN is shown 
in Figure [^. A corresponding Raman spectrum is shown 
in Figure^^ revealing the typical 2D line shape, which 
consists of four Lorentzians (see inset in Figure 

For the low temperature magneto-Raman measure¬ 
ments we employ a commercially available confocal Ra¬ 
man setup, allowing us to perform spatially-resolved ex¬ 
periments at a temperature of 4.2 K and magnetic fields 
of up to 9 T. We use an excitation laser wavelength of 
532 nm with a spot diameter on the sample of around 
500 nm. For detection, we use a CCD spectrometer with 
a grating of 1200 lines/mm. All measurements in this 
work are performed with linear laser polarization. 

In Figure we show a color-encoded two-dimensional 
plot of the Raman intensity as a function of magnetic 
field and Raman shift for a spot on the area of SLG 
encapsulated in hBN. For the hBN and 2D line we ob¬ 
serve only a weak B field dependence (see Supplementary 
Information), in agreement with earlier studie^^. The 
G line, however, shows significant B field dependence, 
which result from the resonant coupling of the G mode 
to inter-band LL transitions whose energies are given by 
Tn = Sn+i + Sm where Sn = VY^/^heBn is the absolute 
value of the energy of the nth Landau level. In particu¬ 
lar, the feature at H « ±3.7 T (see arrows in Figure]^) 
can be attributed to the resonant coupling of the Ti- 
transition (see Figure |^) and the G mode. A close-up 
around the G line is shown in Figure [^, with the arrow 
highlighting the Ti-MPR. In contrast to the magneto- 


Raman spectra taken on SLG encapsulated in hBN, the 
spectra obtained on SLG deposited on Si02 do not show 
any resonant behavior as highlighted in Figure [^, in¬ 
dicating the strong influence of the substrate material 
on the properties of graphene. We therefore focus our 
discussion on the hBN-SLG-hBN data. Interestingly, a 
close inspection of the individual Raman peaks at differ¬ 
ent B fields (see Figure]^) reveals that the G line can be 
well described by a single Lorentzian at all B fields even 
across the resonance. 

In Figures and [^, we show the evolution of the 
measured position ujg and width F^ of the G line with 
B field. We observe three resonances at around 2.1 T, 
3.7 T, and 5.8 T (see e.g. arrows in Figure]^). At these 
values, the energies of the LL transitions T 2 , Ti, and Li 
(see below) match the energy of the G mode phonon at 
zero B field /ico’ph- This is highlighted in Figure!^, where 
we show the dependence of and Li on B field (see also 
vertical dashed lines in Figure]^. The resonant behavior 
of the most prominent feature at 3.7 T can be further 
visualized by plotting the data in a three-dimensional 
representation (Figure and as an Argand diagram 
(Figure]^). When projected on the ojq-^g plane (Fig¬ 
ure]^) the resonance appears as a circle. This behavior 
is a hallmark of a special case of the dynamics of a non- 
Hermitian two level systenP^MSl^ The horizontal width 
of the Argand circle ATg ^ 20 — 25 cm“^ depends on 
the strength of the electron-phonon interaction at the Ti 
level crossing as well as on the charge carrier density rigi, 
which might block LL excitations due to the Pauli prin¬ 
ciple (see Supplementary Information), and the lifetimes 
of the involved states. 

In order to separate the influence of these effects and 
to gain quantitative values, an independent control of 
the charge carrier density is desirable. This is achieved 
by repeating the MPR-experiment on an electrically con¬ 
tacted hBN-graphene sandwich device. Although differ¬ 
ent processing steps, including electron beam lithogra¬ 
phy, reactive io n etc hing, and metal (Gr/Au) evapora¬ 
tion, are needecP^^, the sample quality allows the ob¬ 
servation of MPRs of similar quality and at similar B 
fields as for the unprocessed sample in Figure The 
inset of Figure shows the investigated device. A four- 
terminal back gate characteristic of the graphene resis¬ 
tance is shown as a black trace in Figure and a carrier 
mobility of around /i ~ 4 x 10^ cm^/(Vs) and a width 
of the conductance minimunJ® of ~ 7 x 10^^ cm“^ 
are extracted from these data (see Supplementary Infor¬ 
mation). The high quality of our sample is also seen in 
the quantum Hall measurement (see Figure]^), where 
already at H = 3.7 T well-separated LLs are established, 
which is in good agreement with our observation of the 
Ti-MPR. To pin the charge carrier density close to the 
CNP, we make use of the recently reported photo-doping 
effect in graphene-hBN-heterostructureJ^. When shin¬ 
ing light of sufficient intensity on a gated graphene-hBN 
sandwich device, nitrogen vacancies or carbon impurities 
in the hBN can get charged up to the point where they 
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FIG. 2. Magneto-Raman spectroscopy on SLG on hBN and Si02. (a) Raman spectra plotted as a function of perpendicular 
magnetic field at T = 4.2 K for a spot on the hBN-SLG-hBN sandwich area. The G line shows clear features around B = ±3.7 T 
(see arrows), (b) Schematic illustration of the linear graphene band structure, highlighting the different Landau levels and the 
inter-band Landau level transitions Ti and Li. (c) Glose-up of the G line data shown in panel (a), (d) G line as a function of 
magnetic field for a spot on hBN-graphene on Si02. (e) Individual Raman spectra for different values of the magnetic field on 
the hBN-graphene-hBN area. The Raman spectra are vertically offset for clarity. The G line can be well described by a single 
Lorentzian in all cases. 


completely screen the electric field due to the applied 
back-gate voltage. As a consequence, laser illumination 
of a graphene-hBN sandwich device pins the carrier den¬ 
sity in graphene very close to the CNP, independent of 
the applied gate voltage. This happens on time scales 
much faster than the Raman acquisition time, thanks to 
a laser power of 2 mW at a spot size of 0.25 pm^, which 
corresponds to an intensity that is a factor ^10^ higher 
than that used in Ref. |311 The pinning of the carrier 
density results in the fact that the position and the line 
width of the G line become essentially independent of the 
back-gate voltage, as shown in Figure [^. When turning 
off the laser, the CNP remains pinned to the last value 
of the gate voltage that was applied before. This mem¬ 
ory effect, also called photo-induced dopin^^, is shown 
in Figure [^, where the red and blue traces represent two 
back gate characteristics that are shifted by the photo- 
induced doping effect, but are otherwise nearly unmodi¬ 
fied. Importantly, it has been shown that the width of 
the conductance minimum is independent of the photo¬ 
doping valu^^. This quantity is thus a good estimate 
for the upper limit of the pinned charge-carrier density 
in graphene under Raman measurements. For the inves¬ 
tigated sample we find |nei| < 7 x 10^^ cm“^ (see 

Supplementary Information). 

Having knowledge of the charge carrier density, we can 
now follow AndcP^and Goerbig et al^ to quantitatively 
investigate the low B field MPRs (see gray data points 
in Figures]^ and[^). It has been shown that the depen¬ 
dence of the position ujq and the width F^ of the G line 
as a function of the magnetic field can be well understood 
in terms of the renormalization of the phonon propaga¬ 
tor due to electron-phonon interaction. The phonon self¬ 


energy n(co’) for (nearly) zero doping is given by 

/AT N \ 

n(ci;) = ev^B i\Y,aTn;Uj) + \LY.aLn,Uj)], 

\ n=0 n=l / 

where 

= 2(g - _^_ ^ 

{hujy - {e- ifi7ei,£/2)2 5 - ' 

In the equations above, = 2en represents the 
An = 0 inter-band LL transition energies, while de¬ 
notes the energies of the An = 1 inter-band LL tran¬ 
sitions as defined above (see also Figure [^). The ac¬ 
cording 7 ei,Tn and 7 ei,L„ account for the finite lifetimes 
of the different LL excitations and A and \l denote the 
coupling constants of the T^- and Ln-transitions to the 
G mode, respectively. The coupling of the Ln-transitions 
to the G mo de, which has also been observed in previous 
studie^^^^, is only due to higher order processeP^, gen¬ 
erally making Xl X. The Li transition, in particular, 
has to be included as it is needed to account for the res¬ 
onance at around 5.8 T. Since the charge carrier density 
in our sample is low enough not to affect the strength 
of the resonances, we do not need to take into account 
the filling factor dependence of n(co’). Furthermore, near 
zero doping, all intra-band LL excitations are blocked by 
the Pauli principle. 

The above expression for the phonon self-energy can 
finally be used to find the pole of the renormalized 
phonon propagator by solving the equation — 

2ujpiJl{uj) = 0. Its roots are related to the position ujg 
and FWHM F^ of the G mode via uj = ujq — irG'/2. To 
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FIG. 3. Magneto-phonon resonances on hBN-SLG-hBN 
heterostructure, (a) and (b) G line position (a) and G line 
FWHM (b) as functions of B field as obtained from Lorentzian 
fits to the data. The arrows highlight the three visible 
MPRs. (c) Inter band Landau level transition energies Tn 
(blue lines) and Li (red line) as a function of magnetic field. 
The black line represents the G mode frequency at zero mag¬ 
netic field, cjph- For the plot, effective Fermi velocities of 
t'F = 1.17 X 10® m/s (blue curves) and r'p = 1.12 x 10® m/s 
(red curves) and Ct;ph = 1587 cm“^ are used. 




Vq (cm-') 

FIG. 4. (a) Three-dimensional representation of the data 

shown in Figure (b) ujg versus Fg. The B field value 
of each point is color encoded and corresponds to the points 
shown in panel (a). 


account for non-electronic broadening effects, we make 
the replacement cjph ^ cjph — '<^7ph/2, with 7ph being the 
broadening of the phonon due to non-electronic processes 
and cjph the G mode frequency at zero B field. 

In Figures |^, and ^ we show the comparison be¬ 
tween calculation and experiment. The strength of the 
MPRs is comparable to that on the unprocessed sam¬ 
ples, showing that our etching and contacting techniques 


do not reduce the device quality. The red traces show 
the theoretical result taking a constant Fermi velocity 
of Up = 1.17 X 10^ m/s and a constant broadening pa¬ 
rameter 7ei = 7ei,T^ = 7ei,Ln = 1^0 cm“^. These two 
parameters were chosen such that the Ti-MPR is de¬ 
scribed well. Values for the other parameters are ob¬ 
tained as follows. While the G mode phonon frequency 
at 5 = 0 T, cjph = 1586 cm“^, can be directly extracted 
from the data, the non-electronic phonon broadening 
7ph = 5.5 cm“^ can be extracted from the residual F^ 
at higher magnetic fields (e.g. at 8 T) where no LL tran¬ 
sitions are energetically matched with the phonon mode. 
The electron-phonon coupling parameter A = 4 x 10“^ is 
fixed by the value of F^ at very low magnetic fields as 
it has to guarantee the right amount of Landau damp¬ 
ing. The extracted A i s in good agre ement with values re¬ 
ported by other groupFor the phenomeno¬ 
logically introduced Ap we used a value of Ap = 0.015 A. 
A constant 7ei for all LL excitations evidently overesti¬ 
mates the strength of the MPRs at low magnetic field 
values (see arrows in Figure 5d). Moreover, the positions 
of the resonances do not match between theory and ex¬ 
periment if a constant Fermi velocity is used for all LL 
excitations (see vertical dashed lines in Figures and 
1^). In particular, the Li-MPR appears in the calcula¬ 
tion at a lower magnetic field value than in the exper¬ 
imental data. We thus phenomenologically modify the 
model and assign different widths and effective Fermi ve¬ 
locities (following ShizuygP^ to the three visible MPRs 
(see caption of Figure and Table ^ and end up with the 
blue traces in Figures^^-5f, which are in better quanti¬ 
tative agreement with the experimental data. 

Interestingly, we observe that the extracted effective 
Fermi velocities are significantly higher compared to ear¬ 
lier MPR and infrared measuremen ts on graphene and 
graphene related systemj^l ^^ l ^^ l ^^ l ^^ M^. We attribute the 
high values of vp to many-body effects which arise since 
our samples consist of exfoliated graphene single-layers 
with low doping and low doping fluctuations. So fa r only 
suspended SLG showed similarly high value^^^^. The 
Up values we extracted from the measurement of MPRs 
on 13 spots on three different samples are displayed in 
Table The errors on up are below 0.01 x 10^ m/s. The 
extracted values indicate that up decreases with B field 
and might also depend on the transition as has been pre¬ 
dicted by theor3l^3Sni. 

Finally, we find that the widths of the LL excitations 
are strongly dependent on the LL indices of the involved 
states (see last column in Table 1). The electronic life¬ 
time decreases with increasing LL index, resulting in 
strong suppression of all for n > 3. This observa¬ 
tion might be attributed to the additional decay chan¬ 
nels arising from an increased number of energetically 
lower and unoccupied LLs for excitations with higher n. 
However, we cannot gain independent insights on the B 
field dependence of the LL broadenin^^ as we are limited 
to extract 7e/ for each transition from the corresponding 
MPR. 
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FIG. 5. Transport and magneto-Raman measurements of a contacted hBN-graphene-hBN heterostructure, (a) Longitudinal 
sample resistance as a function of back gate voltage after photo doping at I 4 = — 20 V (red), -10 V (blue) and 0 V (black). 
The inset shows an optical image of the device, with the scale bar representing 5 pm. (b) Longitudinal conductivity axx (red) 
and Hall conductivity axy (blue) plotted as a function of charge carrier density at B = 3.7 T. The typical half-integer quantum 
Hall plateaus for graphene are observed, (c) G line position (blue) and G line FWHM (green) as a function of Vg. (d) and 
(e) (jJg and Fg as functions of magnetic field recorded with Vg set to 0 V. The data points are extracted from Lorentzian fits 
similar to Figure 3. A running average of 25 mT was used. The red lines represent a theoretical model which uses constant 
= 1.17 X 10® m/s and constant 7 ei = 160 cm ^ for all resonances. The blue lines represent the theoretical model, where 
different effective Fermi velocities {vf,Ti = 1.17 m/s, i’f,T 2 = 1.19 m/s, vf,Li = 1.12 m/s) and electronic broadening parameters 
( 7 ei,Ti = 160 cm“^, 7 ei,T 2 = 270 cm“^, 7 ei,Li = 80 cm“^) are employed for the Ti-, T 2 - and Li-resonances. The Tn for n > 3 
are modeled with 7 ei,T^ = 400 cm“^ to fully suppress them, while 7 ei,To was set to 20 cm“^. The other parameters were fixed 
to ct;ph = 1586 cm“^, 7 ph = 5.5 cm“^, A = 4 x 10“®, and Al = 0.015A. The inset of panel (e) and panel (f) show close-ups of 
the T 2 - and Li-MPRs, respectively, highlighting the differences between the blue and red traces. 


Sample: 

#A 

#B 

I’F (xlO® m/s) 

#C 

#C 

7ei (cm 

T 2 (~ 2.1 T) 

1.20 

1.19 

1.20 

270 

Ti 7 3.7 T) 

1.17 

1.16 

1.17 

160 

Li 7 5.8 T) 

1.13 

1.13 

1.12 

80 


TABLE 1. Effective Eermi velocities extracted from the position of the MPRs on three different samples. #A: sample shown 
in Eigure^. #B: second, non-contacted hBN-graphene-hBN sample. contacted hBN-graphene sandwich characterized 

in Eigurel^ The values in the first three columns contain the mean values of the effective Eermi velocities as measured on a 
varying number of spots per sample. The last column contains the LL excitation widths used in the calculation. 


Our fabrication technique, yielding samples with high 
carrier mobility and low doping, can also be applied to 
bilayer graphene. This allows us to investigate MPRs 
in exfoliated BLG, which could not have been studied 
earlier. The most relevant LL transitions in BLG are 
inter-band excitations in the lower subbands, which we 
also denote by = ^n+i where Sn now denotes the 


energy of a Landau level in the lower subband of BLG 
(see Supplementary Information) and n > 1. Their evo¬ 
lution with B field is shown in Figure l§i. As for the 
case of SLG, we fit single Lorentzians to every Raman 
spectrum. The extracted values of ujg and F^ as a func¬ 
tion of B field are shown in Figures and [^, respec¬ 
tively. Two resonances at around 6.4 T and 8.1 T can be 
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FIG. 6. Magneto-phonon resonances in bilayer graphene 
encapsulated in hBN. (a) inter-band Landau level transition 
energies Tn (blue lines) as a function of B field. The black 
line represents the phonon frequency at zero B field, (b) and 
(c) ojg and Fg as functions of B field. The data points are 
extracted from Lorentzian fits to data taken on the sample 
shown in Figure [^. A running average of 25 mT was used. 
The red lines are computed using equation Eq. S2 (see Supple¬ 
mentary Information) with the parameters given in the text. 


well identified, which can be attributed to the LL transi¬ 
tions T 5 and T 4 , respectively. The red curves are calcu¬ 
lated using Eq. S 2 (see Supplementary Information), with 
= 1.07 X 10^ m/s, A = 3.5 X 10-^ 7ph = 6.5 cm-^, 
<^ph = 1585.9 cm“^ and 7 ei = 270 cm“^. Overall, the 
magnitudes of the resonances below 9 T are significantly 
lower as compared to the case of single-layer graphene 
(see Figure]^, which, however, is in good agreement with 
our calculation. Interestingly, in contrast to single-layer 
graphene, a constant Fermi velocity is sufficient to de¬ 


scribe our data. Moreover, compared to SLG, the Fermi 
velocity is reduced, indicating that electron-electron in¬ 
teraction effects play a minor role in bilayer graphene. 
However, as we do not know the charge carrier doping in 
our bilayer sample, a suppression of many-body effects 
might also be due to finite charge carrier density. It is 
possible to describe both resonances with the same LL 
exciton broadening. 

In summary, we investigated the B field dependence of 
the G line phonon renormalization of exfoliated single¬ 
layer and bilayer graphene for low magnetic fields. For 
single-layer graphene we have compared the substrate in¬ 
fluence of Si 02 and hBN. While no distinct features could 
be observed on Si 02 , the Raman spectra on hBN showed 
distinct magneto-phonon resonances. For a quantitative 
discussion of the parameters involved in the coupling of 
Landau level transitions to the Raman G mode we in¬ 
vestigated an electrically contacted graphene sample to 
pin the charge carrier density. We found high effective 
Fermi velocities of up to 1.20 x 10^ m/s, which varied for 
the three observed LL excitations. Due to the low charge 
carrier density, we attribute this finding to many-body 
contributions to the energy of the LL excitons. We show 
that Landau level excitations with higher Landau level 
index have reduced lifetimes, possibly linked to the exis¬ 
tence of an increased number of intermediate LL states 
that lead to additional decay channels. Finally, we ob¬ 
served MPRs in exfoliated bilayer graphene. Here, a con¬ 
stant Fermi velocity and a single Landau level excitation 
width for all resonances is sufficient to reach good agree¬ 
ment with theory. Our work paves the way towards a 
more refined understanding of electronic many-body ef¬ 
fects in the presence of magnetic fields in single-layer and 
bilayer graphene and can be easily extended to other two- 
dimensional materials. 
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